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We discuss violation of Bell inequalities by the regularized Einstein-Podolsky- 
Rosen (EPR) state, which can be produced in a quantum optical parametric down- 
conversion process. We propose an experimental photodetection scheme to probe 
nonlocal quantum correlations exhibited by this state. Furthermore, we show that 
the correlation functions measured in two versions of the experiment are given 
directly by the Wigner function and the Q function of the EPR state. Thus, the 
measurement of these two quasidistribution functions yields a novel scheme for 
testing quantum nonlocality. 

PACS: 03.65.Bz, 42.50.Dv 

1 Introduction 

The work of Einstein, Podolsky, and Rosen (EPR) has started a lasting debate about 
the completeness and the meaning of local realities in quantum mechanics. In order to 
make their point, EPR used the following wave function for a system composed of two 
particles [1]: 



Compared to the original paper of EPR, we have put here spatial separation between 
the particles equal qo = 0. The above wave function describes an entangled state of two 
particles a and b, with the following properties: 

[qa - (Zh)|EPR) = 0, (pa + Pb)|EPR) = 0. (2) 
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However, further studies of quantum nonlocality and entanglement, especially those 
providing quantitative tests of compatibility of quantum mechanics with local theories 
in the form of Bell inequalities, used mainly spin-1/2 particles instead of systems with 
continuous degrees of freedom. 

Quantum correlations for position-momentum variables associated with the EPR 
state (1) can be analyzed in phase space using the Wigner function or the positive Q 
function. Using the Wigner function approach. Bell [2] has argued that the EPR wave 
function (1) will not exhibit nonlocal effects, because its Wigner function is positive 
everywhere, and as such will allow for a hidden variable description of the correlations. 
This statement is correct as long as the measured observables have a straightforward 
classical interpretation in the phase space representation. However, the situation can 
change dramatically, if we take into account quantum observables for which the phase 
space cannot serve as a local theory model. In a recent piiblication [3] we have demon- 
strated that the Wigner function of the EPR state, though positive definite, provides 
a direct evidence of the nonlocal character of this state. The demonstration was based 
on the fact that the Wigner function can be interpreted as a correlation function for 
the joint measurement of the parity operator. 

In this presentation we review the EPR state in the Wigner representation and ex- 
tend our approach to the positive Q function, which also can be regarded as a correlation 
function for the joint measurement of certain dichotomic observables. The analysis will 
be performed for an optical realization of the EPR state as the entangled two-mode 
state of light generated in the spontaneous parametric process. In this quantum op- 
tical context, analysis of the Q function is of particular interest, as the experimental 
demonstration of nonlocal correlations exhibited by the Q function poses less stringent 
requirements on single photon detectors used in the setup. 

This paper is organized as follows. First, in Sec. 2, we briefly review the quantum 
optical version of the EPR state. In Sec. 3 we discuss the Wigner and the Q functions 
of this state, and study the limiting case in which the original EPR state is recovered. 
In Sec. 4 we present the quantum optical scheme which can be used to demonstrate 
nonlocality of the EPR state in the phase space. Finally, Sec. 5 concludes this presen- 
tation. 

2 EPR state in quantum optics 

Before we show that the phase space representation of the EPR state fully exhibits 
the nonlocal character of this entangled state, we shall give a brief description of an 
optical realization of this state, in terms of the two-mode correlated light generated in 
the spontaneous parametric down-conversion process. 

With a clear application to quantum optics, it has been shown that a state produced 
in a process of nondegcnerate optical parametric amplification (NOPA) [4] , is the op- 
tical analog of the EPR state in the limit of strong squeezing. The NOPA process is 
a nonlinear interaction of two quantized modes (denoted by a and b) in a nonlinear 
medium with a strong classical pump field. The strength of the interaction can be 
characterized with the parameter x, which involves the second-order susceptibility and 
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the pump field amplitude. Tlie interaction Hamiltonian of the system is: 

H = ixia^'b'' - ah). (3) 

If the initial state of the system consists of two vacuum modes the NOPA generates: 
|NOPA) = e''(^^^^-'*^)|0,0), (4) 

where r = is a dimensionless parameter characterizing the interaction time. Simple 
algebra, based on the following disentanglement of the evolution operator: 

^{a^y^ — ab) gtanhra^S^ | | tanhrdS 

\ cosh r / ' 

shows that the generated state has a diagonal decomposition in terms of the number 
states of the two modes: 

oo 

|NOPA) = — J— ^(tanhr)"|n,n). (6) 

In order to sec the relation between this state and the EPR state (1), we rewrite the 
state (6) in the following form: 

INOPA) = f;(tanhr)" j dq j dq' \q,q'){q,q'\n,n). (7) 

Using the fact that the scalar products can be expressed in terms of the Hermite poly- 
nomials (in dimensionless units): {q\n) = (2"n!0r)~^/^-ff„(g) exp(— and the 
following summation formula (valid for arbitrary A < 1): 

we obtain: 

, ^ .V 1 /■ f , ( + - 2(7(7' tanhr\ , , , 

This formula is a regularized version of the EPR state (1), with a gaussian smoothing 
profile of the plane waves. Now it becomes clear, that in the limit of r — > oo i.e., for a 
very long interaction time, this smoothing fimction becomes a sharp function: 5{q — q'). 
This follows from the fact that for A = 1, the sum (8) reduces to the completeness 
relation for the oscillator wave functions. In this limit the state (9) becomes exactly 
the EPR state (1), and as result we obtain, that in terms of photons, the EPR state is: 

|EPR) ~ lim INOPA) ~ |0, 0) + |1, 1) + |2, 2) + . . . . (10) 

r— ^oo 

The NOPA state has been generated experimentally [5] and applied to discuss the 
implications of the positivity of the corresponding Wigner function on the Bell inequality 
[6]. 
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3 EPR state in phase space 

In this section, we discuss the Wigner and the Q functions of the state produced in the 
nondegenerate spontaneous parametric down- conversion process. We shall pay partic- 
ular attention to the limit of strong interaction for r — > oo, where the original, singular 
EPR state is recovered. 

3.1 The Wigner function in the coherent state representation 

The two mode Wigner function of the NOPA state (9) can be calculated directly from 
the definition: 

r (fa' r (PB' 

VF(a;/3) = J — J ^exp(aa'* - a*a' + f3p'* - (3'){D^{a')D,{p')) , (11) 

where Da{ct) and Db{P) are the Glauber's displacement operators for modes a and b. 
Simple calculation shows that the Wigner function of the state (9) has the form: 

W{a;P) = ^exp[-2cosh2r(|a|2 + |/3|2) + 2sinh2r(a/3 + a*/3*)]. (12) 

TT 

This positive everywhere Wigner function is plotted in Fig. 1 for real values of a and 
P and for r = 1. The Wigner function of the original EPR state (1) is obtained in the 
limit r — > oo: 

W{a;P) ^ d{ar- f3r)Sia,+ f3,), (13) 

where ar (a^) and Pr {Pi) are the real (imaginary) parts of a and p. Note the singular 
form of this Wigner function, due to the singular character of the original EPR wave 
function (1). 

3.2 The Q function in the coherent state representation 

The two mode Q function of the NOPA state can be calculated directly from its defini- 
tion: 

Q{a;p) = ^\{a,p\NOPA)f. (14) 
It is easy to show that for the NOPA state (9) this function has the following form: 

= 7 Kr^ exp (-lap - + tanh r(a*/3* + a/3)) , (15) 

(TTCOshr)^ 



with the following marginals: 
1 

7r(coshr)^ ^"''^ \^ (coshr) 

1 ^^^p ( \P\' 
7r(coshr)2 \ (coshr) 



Q{a) = — — ^ „ exp 

~ y ^^^;osJ 

= exp (-r^^ ] . (16) 
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Fig. 1. Plot of the Wigner function of the regularized EPR state with r = 1 for real values of 
a and (3. 

This two-mode Q function is plotted in Fig. 2 for real values of a and /3 and for r = 1. 
The Q function of the original EPR state (1) obtained in the limit r — > oo is proportional 
to: 



The last term in the exponent is responsible for the correlation of the entangled EPR 
state. Note that if we consider the Q function given by Eq. (17) as a limit of Eq. (15), 
there is a vanishing prefactor with the magnitude proportional to e"^*". The function 
(17) can be derived directly from the EPR wave function (1): 



This simple integration of the coherent state wave functions, in position representation, 
can be calculated and as a result we reproduce the result (17). 



In this presentation we propose to probe the correlations of the NOPA state with a 
scheme involving two photon counting detectors for modes a and b. The setup to demon- 



Q{a; 13) ~ exp(-|ap - |/3p + a*j3* + a(3). 



(17) 




(18) 



4 Testing quantum nonlocality in phase space 
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Fig. 2. Plot of the Q function of the regularized EPR state for real values of a and /3, for 
r = 1. 

strate quantum nonlocality for the Wigner function and the Q function is presented in 
Fig. 3. The photons from the NOPA source are measured by two photon counting 
detectors preceded by two beam spUtters. These beam spUtters, with the transmission 
coefficients very close to one, mix the NOPA photons with two highly excited coherent 
states. It has been shown elsewhere [7], that in the limit of the transmission tending 
to one, the effect of the beam splitters is equivalent to phase space shifts for the two 
modes by -Da (a) and Db{f3), where a and (3 are complex amplitudes of the reflected 
coherent fields. A single-mode version of such an experiment, involving measurements 
of the displaced photon number statistics for simple classical states of light, is reported 
in these proceedings in the context of quantum state measurement [8] . 

We shall show that if the detectors placed in the two arms of the proposed mea- 
surement scheme can resolve the number of absorbed photons, the phase space Wigner 
function can be determined directly. Moreover, it will become apparent that the Wigner 
function measured in this setup describes correlations between the parity of the number 
of photons registered by the two detectors. As the parity is a dichotomic ±1 variable, 
the Wigner function can be therefore directly inserted into appropriate Bell inequalities 
in order to test the nonlocal character of the NOPA state. 

The most efficient single-photon detectors available currently are avalanche photo- 
diodes operated in the Geiger regime. These detectors cannot resolve the number of 
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Fig. 3. The experimental scheme proposed to test nonlocality exhibited by the quantum 
optical realisation of the EPR state. The nonlinear crystal NL pumped by the coherent field 
generates two correlated beams. The high-transmission beam splitters with auxiliary coherent 
fields placed in the paths of these beams perform the displacement transformations Da (a) and 
DbiP). The phase space displaced beams are monitored by the photon counting detectors PCI 
and PC2. 

simultaneously absorbed photons, and provide only a binary answer saying whether any 
photons have been registered or not [9] . We shall show that use of these detectors in 
our scheme leads in a natural way to the measurement of the two-mode Q function, and 
that the binary answer of the detectors makes the Q function a nonlocal correlation 
function which violates the corresponding Bell inequality. 

4.1 Violation of local reality by the Wigner function 

In order to relate the Wigner function to photocount correlations, we will consider the 

case when the detectors are capable of resolving the number of absorbed photons. Let 
us assign to each event or —1, depending on whether an even or an odd number of 
photons has been registered. In this case the joint correlation for the photon number 
parities is: 

E{a; /?) = (NOPA|n„(a) n6(/3)|NOPA), (19) 
where 

tia{a) = Daia){-lf^Di{a), fltW) = (20) 

are parity operators of the modes a and b characterized by the corresponding photon 
number operators fia and fii,. Equivalcntly, E{a; (3) can be written as the overall parity 
operator displaced in the two- mode Hilbert space by the operation I)J(a)(g)i)|(/?). This 
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observable is known in the phase space representation to provide, up to the normaliza- 
tion constant, the Wigner function of the quantum state [7]. On the other hand, we 
can establish a clear analogy of this scheme with joint spin-1/2 measurements. The 
counterpart of the polarizer orientation is now the coherent displacement, which can be 
set freely in each of the two spatially separated apparatuses. Furthermore, classifica- 
tion of the registered number of photons according to the parity provides a dichotomic 
±1 outcome which is analogous to the spin direction. As a result all types of Bell's 
inequalities derived for spin systems can be used to test the nonlocality of the NOPA 
Wigner function. In Ref. [3] we have shown that the Wigner function (12) leads to a 
violation of the the Bell inequality. Here we will just quote the result. For local hidden 
variable theories the following combination: 

B = E{a'; /?') + E{a'; /?) + E{a; /?') - E{a; (3) (21) 

should satisfy the Bell inequality —2<B<2. We have found that for a certain 
selection of coherent displacements the value of this combination in the limit r ^ oo is: 
B = 1 + 3- 2~*/^ Ri 2.19, which clearly violates the Bell inequality. 



4.2 Violation of local reality by the Q function 

Based on the previous result we shall demonstrate how nonlocality of the EPR state 
can be revealed using the positive definite Q quasidistribution function. The quantum 
optical experiment we shall propose does not require detectors that can resolve the 
number of photons triggering the output signal. We will be interested in events when 
no photons were registered. This measurement is described by the following observable: 

da = limje)"" = :e-"»: = |0a)(0„|, (22) 

e— >0 

for the mode a, and a similar observable Ot for the mode b. As we discussed in the 
previous subsection, these observables can be shifted in phase space using auxiliary 
coherent fields and high- transmission beam splitters. As a result of this shifting we 
obtain: 

da{a) = Da{a)6aDl{a), 6^ = Db{f3)dbDl{p). (23) 

Obviously, these observables describe projections on a coherent state |a) (a] for the mode 
a and |/3)(/?| for the mode b. Therefore, statistics of no-count events yields directly the 
Q function of the measured field. 

The quantum mechanical probability of no-count events in both the detectors is: 

Pab{a;p) = (NOPA|d„(a)®66(/3)|NOPA) 

= j^^^^exp{-\a\^ -\(3\^ +tanhr{a*p* +ap)) (24) 

where a and /3 are the coherent displacements for the modes a and b, respectively. The 
probabilities of no-count events on single detectors are: 

Ma) = (N0PA|Q„(a)^4|N0PA) = ^exp(-J^), 
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Fig. 4. Plot of the combination CTl as a function of the displacement intensity J and the 
interaction parameter r. Only values exceeding the upper bound imposed on CTi. by local 
theories are shown. 



(NOPA|7„0d6(/3)|NOPA) 



1 



(cosh r)' 



■ exp 



(cosh ry 



(25) 



We recognize that these probabilities are equal up to a normaUzation factor the two- 
mode Q function (15) and its two marginals for the NOPA state. 

The correlation function Pab{oi,(3) describes a binary or 1 measurement on the 
modes a and b with adjustable parameters of the apparatuses represented by a and 
/?. As a result of this, the Bell inequality derived for a measurement of local realities 
bounded by and 1, can be applied to test the nonlocal character of the NOPA state 
(6), using the Q function. We shall consider the the Clauser-Horne combination [10] 
for a selected set of four displacements: 



Cn = pab{0; 0) + pab{a; 0) + Pab{0; /?) - Pab{a; (3) - Pa(0) - Pb{Q), 



(26) 



which for local theories satisfies the inequality — 1 < CTi < 0. We will take the coherent 
displacements to be real with a = —(3 = \fj ■ For these values we obtain 



CH 



1 

(cosh r)^ 



(2e-^ 



-2J(l+tanhr) _ ^ 



(27) 
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As depicted in Fig. 4, this result violates the upper bound imposed by local theories. 
In the limit of r — > oo and small values of the displacement, this function becomes: 

Cn « SJe-^'', (28) 

i.e., it violates the Bell inecjuality for all values of r, but the violation becomes infinites- 
imally small. This results from the fact that the Q function (17) in this limit is dumped 
by the prefactor 1/ cosh'^ r, decreasing like e~^^. 

5 Conclusions 

In this presentation we have shown that nonlocality of the EPR state can be revealed 
using its phase space representation of the form of the Wigner or Q quasidistribution 
functions. This is possible owing to the observation that these quasidistribution func- 
tions describe nonlocal correlations that can be detected in a certain quantum optical 
scheme. Of course, this scheme is much more general. It can be applied to measure an 
arbitrary two-mode state of light, and the corresponding Wigner and Q functions will 
always have the operational meaning of nonlocal correlation functions [11]. 
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